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Abstract 

Suzaku observation of the edge-on spiral galaxy NGC 4631 confirmed its X-ray halo extending out to 
about 10 kpc from the galactic disk. The XIS spectra yielded the temperature and metal abundance for 
the disk and the halo regions. The observed abundance pattern for O, Ne, Mg, Si and Fe is consistent 
with the metal yield from type II supernovae, with an O mass of about 10 6 Mq contained in the halo. 
These features imply that metal-rich gas produced by type II supernova is brought into the halo region 
very effectively, most likely through a galactic wind. Temperature and metal abundance may be affected 
by charge exchange and dust. An upper limit for the hard X-ray flux was obtained, corresponding to a 
magnetic field higher than 0.5/xG. 
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1. Introduction 

For the understanding of the chemical evolution of 
galaxies and clusters of galaxies, precise knowledge about 
the metal production from Type la (SN la) and Type 
II (SN II) supernovae (SNc) is of vital importance. X- 
ray imaging spectroscopy of supernova remnants (SNRs), 
galaxies and intra-cluster medium (ICM) has shown that 
metal abundances generally vary from source to source. 
However, Sato et al. (2007a) determined the abundances 
of O, Ne, Mg, Si, and Fe in several galaxy clusters based on 
Suzaku observations, and found that the abundance pat- 
terns are commonly represented by a combination of type 
la and II supernovae products with an occurrence num- 
ber ratio of 1 : 3.5, based on the theoretical yields from SN 
II and SN la. In order to further look into the past his- 
tory of the two types of supernova, we need more precise 
knowledge about the metal yields from the different SNe 
types. It is, however, quite difficult to extract pure super- 
nova products, in particular for the X-ray emitting hot 
gas. When we observe young supernova remnants (SNR), 
there is always a mixture of supernova ejecta and sur- 
rounding ISM and line intensities also strongly depend on 
the ionization condition. One way to provide a good con- 
straint is to observe X-ray halos around starburst galaxies, 
which are considered to be maintained by the enhanced 
SNe II activity in the recent time. 

Recent X-ray observations of starforming galaxies 
showed that metals are contained in the extended hot halo 
of the galaxies. For M82 and NGC 253, RGS spectra for 
several sliced regions along the outflow axis showed lines 
from highly ionized O, Ne, Mg, Si and Fe (Read & Stevens 
2002; Bauer et al. 2007). The observed intensity ratios of 



the lines indicates that the gas is cooling as it travels out- 
ward from the galaxy disk and that the gas around NGC 
253 could be partly out of ionization equilibrium. Suzaku 
observations of a "cap" region of M82, which is 11.6 kpc 
north of the galaxy and is possibly a termination region of 
the hot-gas outflow, showed a spectrum consisting of emis- 
sion lines from O through Fe (Tsuru et al. 2007). These 
spectral features strongly suggest that fresh metal-rich gas 
produced in the starforming region is flowing out mainly 
along the minor axis of galaxies. However, metal abun- 
dances in the halo gas has so far not been well-constrained. 
RGS data are limited in statistics, and both the EPIC and 
the ACIS instruments do not have sufficient energy reso- 
lution below 1 keV (e.g. Tullmann et al. 2006). Suzaku 
offers a good opportunity for measuring the mctallicity of 
the outflowing hot gas. 

NGC 4631 is a nearby Sc/SBd galaxy with an edge-on 
morphology with the distance estimated to be about 7.5 
Mpc, where V corresponds to 2.2 kpc. Estimated mass 
by the Tully-Fisher relation is 2.6 x 10 10 M o (Strickland 
et al. 2004). The inclination and position angle are 81° 
and 356° respectively This galaxy is suitable for Suzaku 
observations of the X-ray halo maintained by its SN activ- 
ity. With its radio halo (Hummel & Dettmar 1990) and 
warm IR ratio, it is classified as a mild disk-wide star- 
burst galaxy (Golla & Wielebinski 1994). An extended 
X-ray halo was discovered by ROSAT (Wang et al. 1995), 
and it has been well studied with Chandra (Wang et al. 
2001; Oshima 2003; Strickland et al. 2004) and XMM- 
Newton (Tullmann et al. 2006). The size of the halo is 
several arcmin and no X-ray counterpart for the central 
AGN has been detected (Strickland et al. 2004). The as- 
sociation of an Ha filament with the X-ray emission has 
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been discovered (Wang et al. 2001), and FUSE also de- 
tected O VI lines from a region at 2' above the disk (Otte 
et al. 2003). These strongly suggest that a halo around 
NGC 4631 is the site of galactic outflow or fountain, where 
the gas is floating up from the disk by the SN energy input 
and possibly cooling. 

Another important feature of the halo is the extended 
synchrotron radio emission, which is observed from several 
star-forming galaxies (e.g. Veilleux et al. 2005). These ra- 
dio halos are populated with relativistic electrons together 
with magnetic fields of an order of 10/iG. The radial ori- 
entation of the magnetic field lines suggests that the field 
has been carried into the halo region by the outflow of hot 
gas (Golla & Hummel 1994). In some cases, the hot gas 
may be confined in the halo when the magnetic pressure 
exceeds the thermal pressure (Wang et al. 1995). To con- 
strain how the non-thermal energy is distributed in the 
halo, it is important to look into the possibility of non- 
thermal emission in the hard X-ray band. 

Throughout this paper we adopt the Galactic hydro- 
gen column density of Ah = 1-27 x 10 20 cm' 2 (Dickey &; 
Lockman 1990) in the direction of NGC 4631. The solar 
abundance table is given by Anders & Grevesse (1989), 
and the errors are the 90% confidence limits for a single 
interesting parameter. 

2. Observation and Data Reduction 

2.1. Observation 

Suzaku carried out observations of the starburst galaxy 
NGC 4631 in November 2006 (PI: N. Y. Yamasaki) with 
an exposure of 81 ks. The observation log is given in 
table 1. The X-ray Imaging Spectrometer (XIS: Koyama 
et al. 2007) instrument consists of a back-illuminated (BI: 
XIS 1) CCD sensor and two front-illuminated (FI: XIS 0, 
3) sensors. The XIS was operated in the Normal clocking 
mode (8 s exposure per frame), with the standard 5x5 
or 3 x 3 editing mode. The observed X-ray image in 0.6 - 
0.7 keV band and XIS contour image in 0.5-2.0 kcV range 
overlaid on the optical image by DSS is shown in figure 1 

2.2. Data Reduction 

We used version 2.0.6 processing data (Mitsuda et al. 
2007) , and the analysis was performed with HEAsoft ver- 
sion 6.4 and XSPEC 11.3.2aj. In the analysis of XIS data, 
we selected ELEVATION > 15° of the data set to remove 
stray-light from the day Earth limb, and the light curve of 
each sensor in the 0.3-10 keV range with a 16 s time bin 
was also examined so as to exclude periods of an anoma- 
lous event rate greater or less than ±3a around the mean. 
After the above screenings, the remaining exposures of the 
observations decreased by ~ 7% as shown in table 1. Event 
screening with cut-off rigidity (COR) was not performed 
in our data. 

In order to subtract the non X-ray background (NXB), 
we used a dark Earth database, provided by the "xisnxb- 
gen" ftools task (Tawa et al. 2008). Although it is known 
that the optical blocking filters of the XIS have grad- 



ually been contaminated by outgassing from the satel- 
lite, we included these effects in the calculation of the 
Ancillary Response File (ARF) by the "xissimarfgen" 
ftools task (Ishisaki et al. 2007). We then generated 
two different ARFs for the spectrum of each region, A v 
and A B , which respectively assumed uniform sky emis- 
sion and the observed XIS image. The task wrote a value, 
"SOURCE.REG.RATIO", which represents the flux ratio 
in the assumed spatial distribution on the sky inside the 
spectral accumulation region to the entire model, in the 
ARF file header to evaluate the flux within the region. 
Since the energy resolution has slowly degraded after the 
launch, due to radiation damage, this effect was included 
in the Redistribution Matrix File (RMF) by the "xisrmf- 
gen" ftools task. 

We reprocessed the HXD data with the CALDB 
files of 2008-04-01 version, and applied data reduction 
by standard criteria. The HXD data was cleaned by 
the Earth elevation angle > 5°, the cut-off rigidity 
(COR) > 6 GV, and was also processed by apply- 
ing the dead-time correction. After the screening, 
the exposure time of the HXD data was 73.3 ks. 
We used the ae_hxd_pinxinome3_20070914.rsp file 
and ae801019010hxd_pinnxb_cl.evt as the response 
and the non X-ray background (NXB) files, respec- 
tively, and also simulated the 100 Ms accumulation 
of the cosmic X-ray background data as described in 
http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/pinjcxb.html 
using ae_hxd_pinflate3_20070914.rsp as a response file. 

3. Analysis and Results 

3.1. XIS spectra and the Galactic background 

We extracted spectra for three regions from the optical 
and X-ray images as shown in figure 1 : 

1. Disk component: a rectangular region of 12' x 3' and 
r < 6' circle centered on (RA, Dec) = (l h 42 m 08 s , 
+32°32'29") (35.5 arcmin 2 ). 

2. Halo component: a circular region of r < 6' from 
the same center position of the disk component and 
outside the disk region (77.6 arcmin 2 ). 

3. Background: outside of the above disk and halo re- 
gion (r > 6'). 

We first fitted the spectra of the background region 
and produced the extra-galactic cosmic X-ray background 
(CXB) and Galactic emission. This was because the O VII 
and O VIII lines from the Galactic emission affected these 
lines from NGC 4631 (see also Sato et al. 2007a). We 
assumed cither one or two temperature apec model for 
the Galactic emission (Lumb et al. 2002), and tested 
the following two models: apec +phabs x power-law, and 
apec-Y +phabs x (apec 2 + power-law), where the apec mod- 
els had fixed metal abundance at 1 solar with zero redshift 
and the absorption column for "phabs" was fixed to the 
Galactic value. The 0.4-5.0 keV spectra for the BI and 
FI sensors for all regions were fitted simultaneously, ex- 
cluding an energy range of anomalous response around 
the Si K-edge (1.825-1.840 keV). The results of those fits 
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Table 1. Suzaku Observation logs for NGC 4631. 



Object 


Seq. No. 


Obs. 


date 


(RA, Dec)* 


Exp. 


After screening 










J2000 


ks 


(BI/FI) ks 


NGC 4631 


801019010 


2006-11-28 03:23 


- 2006-11-29 21:55 


(12 h 42 m 09?3, +32°32'52") 


81.1 


75.9/76.1 



* Average pointing direction of the XIS, written in the RA_NOM and DEC_NOM keywords of the event FITS files. 




12 h 42 m 40 s 12 h 42 m 20 s 12 h 42 m 12 h 41 m 40 s 12 h 42 m 40 s 12 h 42 m 20 s 12 h 42 m 12 h 41 m 40 s 

R.A.(2000) R.A.(2000) 



Fig. 1. LeftThe X-ray image of NGC 4631 observed with Suzaku in the 0.6 - 0.7 keV energy band. Right: The X-ray contour 
map in linear scale from 0.5—2.0 keV overlaid on an optical image taken by DSS. For both images, the observed XIS0,1,3 images 
were added on the sky coordinate after removing each calibration source region, and smoothed with a = 16 pixel or ~ 17" Gaussian. 
Estimated components of extragalactic X-ray background (CXB) and instrumental background (NXB) were subtracted, and the 
exposure was corrected, though vignetting was not corrected. The region where energy spectra were extracted are shown by solid 
and dotted lines in the right figure, for the halo and disk regions, respectively. Two Vertical lines show the region where we took 
the surface brightness profiles in section 3.2. 



are given in table 2. In both cases, the surface bright- 
ness of the CXB component are consistent with the av- 
eraged value of ~ 10 photons cm~ 2 s _1 sr _1 keV" 1 at 1 
keV (Gendreau et al. 1995) with fluctuations for this sky 
area. Afterwards, we adopted the two temperature model 
based on the F-test of these results. The intensity of the 
OVII line is 4.4 ±0.8 photons cm _2 s _1 sr _1 (this unit is 
hereafter refereed to as the Line Unit or LU), which is 
almost the same as or slightly smaller than the previously 
reported values (see e.g. Sato et al. 2007b; McCammon 
et al. 2002). If we subtract the energy spectrum of the 
background from those of the halo and disk regions, the 
resultant line intensities without correcting the absorption 
effect are 17 ' .h + _\\ LU for O VII and U.OtH LU for O VIII 
in the halo and 5.1 LU for OVII and 5.8±g;| LU for 
O VIII in the disk region. 

In order to take into account both the existence of 
the Galactic component itself and propagation of its sta- 
tistical error, we simultaneously fitted all (BGD, halo, 
disk) regions with the same model: phabs x zphabs x 
(vapec 1Tor2T + zbremss + power-law) + constant x (apec 1 + 
phabs x apec 2 )- In the model, phabs corresponds to our 
Galactic absorption with fixed Ah = 1-27 x 10 20 cm -2 . The 
term (apec 1 + phabs x apec 2 ) represents the Galactic com- 



ponent with a normalization factor to keep a uniform sur- 
face brightness. Based on the previous fits, we fixed the 
temperatures of the Galactic emission, apec 1 and apec 2 , 
to be 0.114 keV and 0.310 keV, respectively, with 1 solar 
abundance and zero redshift. The power-law component 
which represents the CXB was set to be common for all 
three regions. 

For the ISM emission from NGC 4631, zphabs corre- 
sponds to the intrinsic absorption. The hot X-ray halo 
around NGC 4631 is modeled either with a one or two tem- 
perature collisionally ionized plasma with various metal 
abundances given by the vapec model. Note that the abun- 
dances were common for both temperature components, 
and we combined the metal abundances into five groups 
as O, Ne, (Mg & Al), (Si, S, Ar, & Ca), and (Fe & Ni). As 
we found that the continuum emission above 2 keV cannot 
be fitted only with the CXB, we added a zbremss model 
of kT = 10 keV for the low-mass X-ray binary (LMXB) 
component. For the spectrum of the background region, 
zphabs and the normalizations of vapec and zbremss were 
all fixed to be 0. 

The resultant spectra are shown in figure 2 and parame- 
ters are summarized in table 3. For the single temperature 
ISM model, the resultant temperatures are kT~ 0.297 keV 
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Table 2. The best-fit parameters of the apec components for the spectra in the background region of NGC 4631 with one or two 
temperature models (apec) for Galactic emissions, and a power-law model for CXB. 



Parameters 


apeci -f- phabs x power-law 


apec 2 -\- phabs x {apec 2 -\- power- law^) 


kT t (keV) 
Norm i * 


n 145+u-uum 

U.l<±O_ 010 

0.80 ±0.13 


114 +o.oiu 
1.02 ±0.42 


kT 2 (keV) 
Norm2 * 

r 

Norrn^ 


0.93 ±0.05 


n Qin+0.096 
U.*UU_ 061 

0.13 ±0.13 
1 58+ - 07 

0.85 ±0.09 


X 2 /dof 


396/375 


382/373 



* Normalization of the apec components divided by the solid angle, Q 1 ' , assumed in the uniform-sky ARF calculation (20' radius), 
Norm= r nc n H dy/(47r(l + z) 2 _D|)/n D xlO 20 cm 5 arcmin 2 , where D\ is the angular distance to the source. 

* Normalization of the power-law component divided by the solid angle same as the normalization of apec, in units of 10 — 8 f2 u — 1 
photons keV - 1 cm -2 s _1 arcmin -2 at 1 keV. 



for the disk and 0.290 keV for the halo, respectively. The 
emission line profiles for Si and S are difficult to distin- 
guish, and the abundance value is almost consistent with 
0. The fit statistics shown in table 3 support the two 
temperature model. The two temperature model gives 
larger absorption column and larger intrinsic luminosity 
than the single temperature model. The power-law index 
of the CXB component was obtained to be 1.48i"qq3 anc ^ 
the intensity of the LMXB component for the disk cor- 
responds to ~ 3 x 10 39 erg s _1 , which is consistent with 
the Lx — M relation for spiral galaxies (Gilfanov 2004). 
A fraction of the LMXB component in the halo region is 
considered to be due to the extended tail of the telescope 
PSF. We consider that the two temperature model actu- 
ally approximates multi-temperature emission, and it is 
represented by emission peak temperatures for O VII and 
O VIII with lower metal abundance. However, we do not 
further carry on detailed multi-temperature analysis given 
the statistical quality of the data. 

We note possible systematic errors due to the assump- 
tion of the contamination model of the XIS. When we 
change the thickness of the contamination by ±10 %, the 
flux and the temperature of the kT\ component of the 
halo changes about 30%, and 3% respectively. Thus, the 
abundances of metals also change by ~ 10 %. The kT 2 of 
the halo and kT\ and kT 2 of the disk change only by 1%. 

We briefly compare with the previous results for NGC 
4631 here. Wang et al. (2001) reported higher tempera- 
tures (0.18 ±0.02 and 0.61 ±0.1 keV) and lower (0.08 ± 
0.04) metal abundances after subtracting the blank field of 
the similar foreground absorption. Tullmann et al. (2006) 
observed with XMM-Newton and detected halo emission 
up to 9.1 kpc from the disk. They sliced the halo into 8 re- 
gions and fitted the spectra with 2-temperature Raymond- 
Smith plasma model with the cosmic metal abundance. 
The resultant temperatures are kT so f t ~ 0.17 — 0.07 keV 
and kThard ~ 0.28 — 0.2 keV. Considering the different 
plasma code and assumed metal abundance, the agree- 
ment of the results is fairly good. 

3.2. Hardness ratio 

We produced count rate profiles in a rectangular region 
of 7.5 x 42 kpc 2 as shown in figure 1 in two X-ray en- 



ergy bands, 0.4 - 0.85 and 0.85 - 1.5 keV. With Chandra, 
Strickland et al. (2004) also showed surface brightness pro- 
file in 0.3 - 1 keV along the minor axis and found that 
an exponential model was preferred over the Gaussian or 
power-law models. We fitted the profile with an expo- 
nential + constant model as shown in figure 3, in which 
the constant component represents the sum of the X-ray 
and non X-ray backgrounds. The fit well represent the 
data, and can not be rejected within a confidence level of 
95%. The constant levels are consistent within errors to 
the same for the north and south sides. 

The obtained scale heights are summarized in table 4. 
These values are larger than those by Strickland et al. 
(2004) and could be due to the difference of the data region 
used in the fit. The scale height values in table 4 are 
consistent with each other to within 90% errors, but those 
in the hard band tend to be smaller than those in the soft 
band for both sides of the disk. This corresponds to a 
spectral softening, consistent with the feature obtained 
by the spectral fit with XMM-Newton (Tullmann et al. 
2006). 

We calculated the hardness ratio for the exponential 
component after subtracting the constant levels with er- 
rors, and the results are plotted in figure 3. If we as- 
sume a single temperature thermal emission for the disk 
with abundances shown in table 3, temperatures of kT = 
0.2,0.3, and 0.4 keV give HR values of 0.268, 0.610, and 
0.809, respectively. Toward the outer region at > 10 kpc 
from the disk, kT is consistent with > 0.2 keV. The small 
or no decline of the temperature may imply that the gas 
is adiabatically expanding into a vacuum with small me- 
chanical work. 

3.3. Search for the hard X-ray emission 

Since NGC 4631 is accompanied by a strong radio halo 
which is likely to be due to synchrotron emission by rel- 
ativistic electrons, one can expect inverse Compton emis- 
sion in the hard X-ray band. We searched for hard X-ray 
emission with the PIN detector of the HXD over the 34' 
fov. After subtraction of the NXB component estimated 
by the standard modeling method in the current process- 
ing, the energy spectrum was consistent with the expected 
CXB level. Due to the current uncertainty of the NXB 
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Energy (keV) Energy (keV) Energy (keV) 



Fig. 2. The panels show the observed spectra after subtracting the NXB component for all regions of NGC 4631 which are denoted 
in the panels, and they are plotted by black and red crosses for BI and FI, respectively. The spectra not corrected for the integrated 
area. They are fitted with the model: phabs X zphabs X (™j)ec 1Tor2T + zbremss + power-law) + constant X (apec 1 + phabs X apec 2 ) as 
shown by black and red lines for the BI and FI spectra. For simplicity, only the model components for BI spectra are shown. Green 
and blue lines are the ISM component by vapec, cyan and grey are the Galactic background emission by apec 1 and apec 2 , magenta 
and orange are the LMXB and CXB component, respectively. The lower panels show the fit residuals in unit of a. 



Table 3. Summary of the best-fit parameters of one or two vapec components by simultaneous fit of all regions. 



ISM IT 



ISM 2T 



Disk 



Halo 



Disk 



Halo 



-0.009 
-0.016 
+2.4 
-2.2 

2.2 x 1CT 13 



+0.027 

-0.053 
■1.2 



+0.025 
0.016 
+2.3 



-0.065 
-0.015 

+ 13.4 
-4.2 

2.3 x 1CT 14 

n 300+ 016 

U.OUU_q 027 
o o+2.2 
z - z -1.3 

1.4 X 1CT 13 

(mis-" 
i-09±J:it 

j+1.14 
'-0.69 
9 i 7 +0.83 
z - x ' -2.17 
,+0.31 
-0.23 



kT x (keV) 
Norm\ 

fluxj (erg cm -2 s _1 ) 

kT 2 (keV) 

Normt, 

flux^ (erg cm" 2 s" 1 ) 
O (solar) 
Ne (solar) 
Mg, Al (solar) 
Si, S, Ar, Ca (solar) 
Fe, Ni(solar) 
extra absorption (cm -2 ) 
LMXB flux t (erg cm" 2 s" 1 ) 
Hot gas luminosity* (erg s" 1 ) 
X 2 /dof 



0.297 
4.5 



0.97] 



7+0.81 
-0.55 
L92 +1.70 



1.85 
0.45 
1.09 

4.0 
2.0 



3.4 

3.0 x 10 
1.9 x 10 



-1.25 
+3.41 
1.01 
+2.55 
0.45 
+0.75 
-0.64 

X 10 20 



-13 
39 



0.290 

0-9-0.8 
1.3 x 10" 13 



2.18 
2.29 



1 34+ 106 

i -°^-0.91 
+5.34 
1.37 
+ 14.06 
1.79 
K qo+20.57 
°- y °-4.80 

0.85_ ; 63 

0-0±S"o x 1Q2 ° 
8.5 x 10~ 14 



0.9 x 10 39 



1198/1114 



0.209 

8-0_4.5 

1.6 x 10~ 13 
0.480t°;°oi 

1.1 x 10 -13 

n 7 o+0.69 
u - '°-0.38 

i - ou -0.72 
-+0.70 



1.26 
0.60 



-0.65 
+2.65 
0.60 
+0.64 



0.93! . 41 
x 10 20 

-13 



8.6 
2.6 x 10 
3.1 x 10 39 



0.096 
4.5 



0.98] 



0.46] 



6Atf 4 7 x 10 20 
7.8 x 10~ 14 
1.6 x 10 39 



1172/1109 



* Normalization of the vapec component scaled with a factor of source_r.ATIO_r.eg / AREA, which is Norm = 
source_ratio_reg J n crl}1 dV / [An (I + z) 2 D\\ xlO -20 cm -5 arcmin- 2 , where D A is the angular distance to 
the source. 



t Flux within the accumulated region between 0.5 and 2 keV 
* Intrinsic luminosity between 0.5 and 2 keV 
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Table 4. Scale height of the X-ray emission 



Scale height (kpc) 



North side South side 



Hard band 



Soft band 



(0.85 - 1.5 keV) 

(x 2 /dof) 

(0.40 - 0.85 keV) 

(x 2 /dof) 



3.01 ± 0.47 
(66.02/65) 
3.85 ±0.58 
(79.24/65) 



3.33± 0.47 
(81.30/65) 
3.72 ± 0.59 
(87.08/65) 





10 

Distance (kpc) 



kT (keV) 
0.6 

0.5 
0.4 
0.3 
0.2 




l LL 




j >.fJ 



10 

Distance (kpc) 



kT (keV) 



Fig. 3. Count rate profile of hard (0.85 - 1.5 keV) and soft (0.4 - 0.85 keV) X-ray bands without subtracting background and the 
hardness ratio in the north (left) and the south (right) side. The region is shown in figure 1. Dashed lines in the count rate profile 
are drawn to show the averaged rate taken from the background region for the spectral fitting. 



estimation, we scaled the normalization of the NXB by 
—3% (Kokubun et al. 2007), and the spectrum was fitted 
by phabs x power-law model with a fixed photon index to 
be 2.0 in 12-60 keV. 

The resultant upper limit corresponding to la level of 
the flux was 5.5 x 10 -12 erg cm' 2 s _1 in 12-60 keV, with 
the estimated contribution from LMXBs with an energy 
spectrum of kT = 10 keV is 4 x 10~ 13 erg cm~ 2 s _1 . The 
Very Large Array (VLA) observations at 1.49 GHz in- 
dicated a radio halo of NGC 4631 with a brightness of 
about 1.22 Jy (Hummel & Dettmar 1990). If we assume 
an inverse Compton process with 2.7 K photons scattered 
by the same relativistic electrons responsible for the ra- 
dio halo, the strength of the magnetic field is constrained. 
Following the prescription given by Harris & Grindlay 
(1979), we derived the lower limit of the magnetic field 
strength to be B > 0.5 /iG, which is consistent with the 
previous estimation of 5 fiG from the polarized component 
of the radio emission (Hummel et al. 1991). 

4. Discussion 

Suzaku observations of NGC 4631 showed significant 
emission lines from an extended halo region, including 
those from O for the first time. We examined an abun- 
dance pattern for the X-ray halo in the form of the num- 
ber ratio to O as shown in figure 4. To make this plot, 
Z/O ratios are determined by two parameter errors from 
the simultaneous fit of all regions as described in section 
3.1. In this plot, the SN la yields were taken from W7 



model in Iwamoto et al. (1999). For SN II, Nomoto et 
al. (2006) gave an average yields for the Salpeter's IMF 
of stellar masses from 10 to 50 M Q with the progenitor 
mctallicity of Z = 0.02. If the metallicity of progenitors in- 
creases to Z = 1, relative abundances of ejecta increase by 
at most 20% for Fe/O. Also the solar abundance template 
of Anders & Grevesse (1989), and the average abundance 
pattern for 4 clusters and groups (Sato et al. 2007a) are 
plotted together. The pattern in clusters is well fitted by 
a combination of SN la and SN II with a ratio of 1:3.5. 
The metal abundance of NGC 4631 halo disagrees with 
that of SN la and the cluster average, but is consistent 
with the SN II yields. In the disk component, Fe shows 
relatively higher abundance and the abundance pattern is 
consistent with the Solar abundance given by Anders & 
Grevesse (1989). 

For the study of the energetics of the gas, we estimate 
the density, mass, and total energy of the X-ray emitting 
gas. As shown in figure 3, the X-ray profile of the halo is 
very smooth, and there is no boundary recognized between 
the disk and the halo. Hereafter, we will not discriminate 
the disk and the halo component and treat them at the 
same time. We assume a simple exponential model of 
n{z) = no x exp(— z/h) with z indicating the distance from 
the galactic plane, within a radius of 10 kpc for the coronal 
gas. Since the scale height of the surface brightness is 
about 3.5 kpc as shown in Table 4, we can adopt the 
density scale height to be h = 7 kpc. A single temperature 
model gave us a lower limit of the density at the disk as 
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n = 2 x 10 3 cm 3 , and the pressure as ijqT = 7 x 10 3 
cm~ 3 K, which are consistent with the previous values by 
Wang et al. (1995); Wang et al. (2001). With this density, 
the total mass of the X-ray emitting gas is 1.3 x 10 8 M Q . 
and the stored thermal energy is 2 x 10 56 erg. Assuming 
an O abundance of 0.8 solar, the O mass in the hot gas in 
both the disk and the halo is ~ 10 6 M Q . 

Since the cooling time in the disk is about ~ 6 x 10 8 
yr assuming the cooling function of Sutherland & Dopita 
(1993), the required energy input rate is 3 x 10 47 erg yr -1 
and a mass transfer rate of ~ 0.2M Q yr -1 , respectively. If 
one employs a flow time to a radius of 10 kpc, it is about 
5 x 10 7 yr and 10 times higher rates for energy input and 
mass transfer are implied. Wang et al. (1995) estimated 
the maximum mass flow rate from the density multiplied 
by the sound velocity to be 1.4M Q yr -1 . Based on the 
UV observation of O VI line, Otte et al. (2003) estimated 
the flow rate to be 0.48 ~ 2.8M© yr -1 assuming a cooling 
flow model by Edgar & Chevalier (1986). Considering 
the differences in the assumed physical process and con- 
dition, mass flow rates are in good agreement around an 
approximate value of 1-M© yr -1 . 

The most likely source of the energy and the material 
in the coronal gas is SNe. The SFR is estimated by the 
FIR luminosity as 3M Q yr" 1 (Strickland et al. 2004). As 
NGC 4631 has an edge-on morphology, the emission could 
be underestimated due to absorption through the disk. 
Persic et al. (2004) proposed another method to estimate 
the SFR using the X-ray luminosity of high mass X-ray 
binaries (HMXB) in the 2-10 keV band. A typical lumi- 
nosity ratio of HMXB/LMXB of 0.2 gives a SFR of 1.2M Q 
yr -1 . If we assume that all the flux above 2 keV comes 
from HMXBs, it gives an upper limit as 6M Q yr -1 . We 
note that a spectral fit with a power-law component with 
r = 1.2 as suggested by Persic et al. (2004) requires a steep 
r = 1.55 CXB component. The two independent estimates 
give a consistent SFR of ~ 3M Q yr -1 , which is almost the 
same as or a little less than the level in our Galaxy, i.e. 
~ 5M yr -1 (Zinnecker & Yorke 2007). Ass uming the 
Salpeter's IMF integrated over 0.1 to 5OM , a SFR of 
3M yr" 1 gives a SN II rate of 7 x 10~ 3 SN yr" 1 . In ad- 
dition to SN II, SN la of 2.5 x 10" 3 SN yr" 1 is expected 
for the total mass of 2.6 x 10 10 Af o and the SFR of 3M Q 
yr" 1 (Sullivan et al. 2006). 

If SNe are the source of the energy and the mass of the 
halo gas estimated above, one SN needs to supply an en- 
ergy of 3 x 10 49 erg and a mass of ~ 20A/ Q including an 
O mass of 0.2M Q . Therefore, 3% of the typical explosion 
energy of 10 51 erg and 20M Q from the ejecta and ambient 
material have to escape into the halo. These values indi- 
cate that the halo gas is produced very efficiently. Since 
one SN produces 0.14 and 1.8M Q of O for SN la and SN 
II, respectively (Iwamoto et al. 1999), a supply of 0.2M Q 
of O seems plausible. The abundance pattern in the halo 
is, however, well represented by the SN II products, al- 
though the above estimated rate between SN la and SN 
II is ~ 1 : 3. It may imply the selective escape of SN II 
ejecta into the halo gas, which can result if the occurrence 
of SN II is concentrated in the starforming region. In 



this case, superbubbles are formed and the metal-rich hot 
gas escapes along a chimney to the halo space (Norman 
&Ikcuchi 1989). 

We must note two possibilities which may bias the 
abundance and temperature estimate. Lallement (2004) 
pointed out the possible contribution of charge exchange 
(CX) processes between the galactic wind and gas clouds 
in the halo. The CX spectrum is dominated by emission 
lines, and it tends to decrease the apparent temperature 
of the halo, kT\ in table 3. The process was seriously 
evaluated for the M82 halo (Ranalli et al. 2008) and M82 
"cap" region (Tsuru et al. 2007), and they found that al- 
most all of the O VII triplet could be produced by CX. In 
the case of NGC 4631, a dust arch with a mass of a few 
time 10 8 Mq was discovered in the halo, but there seemed 
to be no connection between the X-ray emission and the 
dust arch (Taylor & Wang 2003). Thus, the contribution 
from CX might not be so large, considering the low ve- 
locity of the outflow and the low density of the neutral 
material, which is also suggested by the slow decline of 
the X-ray temperature. We must wait for improved spec- 
troscopic observation to distinguish the physical process 
of the emission. 

Another possibility is the role of dust. Dust of silicates 
(Mg2Si04) and other forms with a mass of ~ 5M© SN -1 
can be formed 300-600 days after SN II explosions (Todini 
& Fcrrara 2001). Since it requires more than 10 7 yr to 
evaporate in the low-density environment with n = 10 -3 
enr 3 (Tielens et al. 1994), a significant amount of metals 
may be held in dust. We hope that future high resolution 
X-ray spectroscopy will be able to show the ionization 
condition of plasmas in the galactic winds more precisely. 

5. Conclusion 

We determined the temperature and metal abundance 
of the X-ray emitting halo gas around NGC 4631. The 
total energy, mass, and metals in the halo can by supplied 
by SNe with the currently estimated SFR, if the outflow 
efficiently carries metal-rich gas from the starforming re- 
gions into the halo. The effect of neutral material and dust 
should be taken into account to understand the plasma 
properties in the halo. 
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Fig. 4. Number ratios of Ne, Mg, Si, S, and Fe to O for disk 
and halo regions. Solid, dotted, dashed, and dot-dashed lines 
correspond to the number ratios of metals to O for abun- 
dance patterns of SNe II yield of Nomoto et al. (2006), solar 
abundance by Anders &; Grevesse (1989), cluster average in 
Sato et al. (2007a), and SNe la yield of Iwamoto et al. (1999), 
respectively. 
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